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Abstract. We present here the first analysis of the XMM 
Newton EPIC-MOS data of the galaxy group around 
NGC4839, which lies at a projected distance to the Coma 
cluster center of 1.6 IiJq Mpc. In our analysis, which in- 
cludes imaging, spectro-imaging and spectroscopy we find 
compelling evidence for the sub group being on its first 
infall onto the Coma cluster. The complex temperature 
structure around NGC 4839 is consistent with simulations 
of galaxies falling into a cluster environment. We see in- 
dications of a bow shock and of ram pressure stripping 
around NGC4839. Furthermore our data reveal a displace- 
ment between NGC4839 and the center of the hot gas in 
the group of about 300 kpc. With a simple approx- 
imation we can explain this displacement by the pressure 
force originating from the infall, which acts much stronger 
on the group gas than on the galaxies. 

Key words: Galaxies: clusters: general - intergalactic 
medium, general - Cosmology: observations - large-scale 
structure of the Universe - X-rays: general 
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1. Introduction 

There is now compelling evidence from studies in different 
wavelength bands that the Coma cluster has substructure. 
The search for substructure in the Coma cluster starts 



more than a decade ago: Fitchett & Webster (1987) find 
a double structure in the cluster center, and Mellier et al. 



(1988) find a second peak in the large scale galaxy distri- 
bution around NGC4839, which was further confirmed in 



a study by Merritt & Tremblay (1994). While these studies 



were based on optical data, Briel et al. (1992) using X-ray 



data coming from the all-sky-survey of the ROSAT satel- 
lite, found substructures coinciding with NGC491 1, and 
again, NGC4839. Another study by White et al. ( |l993| ) 



based on a deeper ROSAT observation revealed even more 
structure. In this pa per w e concentrate on the XMM- 
Newton (Jansen et al. 2001 ) observation of the sub group 
around NGC4839, which has a projected distance from 
the Coma cluster of 1.6 h§Q Mpc. In section 2 we review 
briefly the observation and data treatment of the EPIC 
data (see Turner et al. 2001 ) . This is followed by Sections 
dealing with imaging and spectro-imaging. In section 5 we 
perform a spectroscopic analysis in different regions which 
is followed by discussion and conclusion. This paper is one 
of three papers which present XMM-Newton data on the 
Coma cluster (see also Arnaud et al. 2001a; Briel et al. 
p00l| ). 



Throughout this paper, we assume Hq = 50 km s 
Mpc" 1 h^and q Q = 1/2. At a redshift of z = 0.0232, 1' 
corresponds to 38.9 kpc. 



2. Observations 

The subgroup associated with NGC4839 was observed 
for 36 ksec with XMM-Newton in full frame mode with 
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Fig. 1. The 0.5-2.0 keV image of the EPIC-MOS1 and 
EPIC-MOS2 camera. To correct for vignetting we cor- 
rected each indivi dual p hoton with a weighting function 
(see Arnaud et al. 2001b ) . The image is gauss filtered with 
a = 30". The spacing of the contours is linear with a step 
size of 1.17 x 10~ 3 cts/sec/arcmin 2 . The lowest contour is 
at 4.68 x 10~ 3 cts/sec/arcmin 2 . 




Fig. 2. The 2.0-5.0 keV image of the EPIC-MOS1 and 
EPIC-MOS2 camera. To correct for vignetting we cor- 
rected each indivi dual p hoton with a weighting function 
(see Arnaud et al. 2001b ). The image is gauss filtered with 
cr = 30". The spacing of the contours is linear with a step 
size of 2.34 x 10 -4 cts/sec/arcmin 2 . The lowest contour is 
at 1.87 x 10~ 3 cts/sec/arcmin 2 . 



the medium filter. Due to remaining problems with the 
CTI (Charge Transfer Inefficiency) calibration in the ex- 
tended full frame mode of the PN camera we concentrate 
in this paper on the data from the EPIC-MOS1 and EPIC- 
MOS2 cameras. The background intensity of the instru- 
ments varies considerably with time over several orders of 
magnitude. In order to optimize the signal-to-noise of the 
observation we reject time intervals with high background 
count rates. For this selection we binned the counts in the 
energy band 10-12 keV in time intervals of 100 sec. In this 
energy band we hardly detect any astrophysical sources, 
due to the very small effective area of the mirrors. We re- 
ject time bins with more than 15 counts per 100 sec. After 
this selection we obtain exposure times with 30.6 ksec for 
EPIC-MOS1 and 30.9 ksec for EPIC-MOS2. 

3. Imaging 

Fig.|l| and Fig.|| display the images obtained in the band 
0.5-2.0 keV and 2.0-5.0 keV, respectively. The band 5.0- 
10.0 keV is dominated by the particle background. As al- 



rea dy ob served with R OSAT (Briel et al. |1992| White et 
al. 1993; Dow & White 1995| ) we see a tail-like structure 
around NGC4839 pointing away from the Coma cluster. 
The orientation of NGC4839 in the optical coincides with 
the orientation of the tail. The X-ray emission of NGC4839 
seems somewhat disconnected to the diffuse X-ray emis- 
sion coming from the group in the center of the pointing. 



4. Spectro-imaging 

In order to look for spectral variations we calculated a 
hardness ratio image of the data (see Figj^). To account 
for the background, especially the particle background, we 
used another Coma pointing with the central coordinates 
of RA=13h01min37sec Dec=27dl9m52s (J2000), which is 
situated South-East of the Coma cluster at roughly the 
same offset angle than the sub-group. (46 arcmin in com- 
parison to 51 arcmin for the subgroup pointing). We ori- 
entate this background pointing in such a way that it ac- 
counts for the emission of Coma itself plus stray light in 
the pointing for NGC4839. However, this rotation has no 
measurable effect on the hardness ratio image, the features 
with or without this rotation rest identical. This indicates 
that the background we measure is dominated by parti- 
cles. 

As we are not interested here in the emission of dif- 
ferent point sources but only in spectral differences of 
the extended X-ray emission we use a wavelet algorithm 
(the MVM algorithm Rue & Bijaoui |1997[ ) to filter out 
point-like structures for the hardness ratio images. Our 
hardness ratio image is calculated as hr = (image(2- 
5 keV)-image(0. 5-2.0 keV))/ (image(2-5 keV)+image(0.5- 
2.0 keV)). 

We do not take into account any vignetting effects in 
the hardness ratio image since we know from calibrations 
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Fig. 4. The temperature estimates of the four regions with 
Fig. 3. The background subtracted hardness ratio image different background observations. The horizontal axis: 
in the bands 2.-5.0 keV and 0.5-2.0 keV. The image shows Coma background, the vertical axis: Lockman hole back- 
the data only up to a radius of 12 arcmin. Regions outside ground. Error bars are 90% confidence level, 
this area do not show sufficient statistics to give mean- 
ingful results. The contours are the results from a wavelet 
filtering of an image in the band 0.5-2.0 keV. 



that the vignetting functions are parallel to each other for 
energies below 5 keV. 

The resulting image is shown in Fig.||. We can see that 
there are indeed spectral variations across the field-of- 
view. Around NGC4839 we see that the tail-like structure 
has a different colour (is softer) than the regions North or 
South surrounding it (the tail has a mean hardness ratio of 
(hr) = —0.6). The region, which shows a harder spectrum 
in the South ((hr) = —0.4) seems also to have a tail-like 
structure, which is parallel to the tail seen in the original 
image. The core of NGC4839 shows (hr) = -(0.5 - 0.6) 

The appearance of a region with hard spectrum in the 
South-East or North, at the border of the hardness ratio 
image is not statistically significant. We cannot state with 
certainty that these regions are hotter due to an interac- 
tion of the group with the intracluster medium (herafter 
IGM) of the Coma cluster. These "hot" regions are not 
unlikely to be an artifact which is created by low photon 
statistics together with remaining stray light. 

5. Spectroscopy 

For our spectral analysis we selected different regions in 
the pointing, which show spectral differences in the hard- 
ness ratio image. In our spectral analysis we restrict our- 
selves to the energy range between 0.3-10.0 keV. The sum- 
mary of the spectral fitting is given in Tab.]!]. The position 
and extraction regions are displayed in Tab.^. The spec- 
tra for some of the regions are shown in Figj5[ For the 



spectral fitting we used XSPEC (Arnaud 1996|) a nd the 



M EKAL code (Mewe e t al. | 1985| ; Mewe et al. |1986| Kaas- 
tra 1992, Liedahl et al. 1995). The galactic absorption was 
fixed to the value of the 21cm line (nH = 9 x 10 19 cm -2 - 
Dickey & Lockman 1990] ). For the background spectrum 
determination we used the observations performed on the 
Coma background pointing (see above). In order to take 
into account instrumental variations across the field-of- 
view we selected for the background the same regions on 
the detector than for the source spectra. As the particle 
contribution is somewhat detector dependent, it is bet- 
ter to take the same regions on the detector rather than 
the regions corresponding to the Coma emission (the ori- 
entated regions, like for the hardness ratio image), since 
our background is highly dominated by particles and not 
by cluster emission/stray light. In order to estimate the 
influence of background variations we also used the obser- 
vations of the Lockman hole as background. For EPIC- 
MOS1 we used the Lockman hole observations with the 
thin filter, which is not entirely consistent, since all Coma 
pointings were observed with the medium filter. However, 
this will give us an estimate on the intrinsic variations of 
the background. The comparison of the best fit temper- 
ature estimates are shown in Fig.[|. As one can see, the 
temperatures are very similar. 

To account for vignetting we used the approach de- 



scribed in Arnaud et al. (2001b) 
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Fig. 5. The spectra of region 2 (left) and region 3 (right). The black crosses show EPIC-MOS1 data, the red crosses 
show EPIC-MOS2 data. 



Region description 



AT 



metal 
Abund. 



X^red. 



3.1? 



core of 
NGC4839 

tail of 
NGC 4839 
region south 

of tail 
main body of 4.4|;J 
sub group 



4.81:8 o.7og:g x 



8.0i 



0-368;lg 



0.82 
0.99 
0.69 
0.87 



<°1.69 



1 M 1 - 63 



Table 1. Summary of the spectral fits in different re- 
gions for the combined spectra of EPIC-MOS1 and EPIC- 
MOS2. The errors are 90% confidence level. In order to 
account for background variations we allowed a 10% er- 
ror in the normalization of the background spectrum for 
MOS1 and MOS2. For the fitting of a thermal plasma we 
used the MEKAL code in XSPEC. 



5.1. Region 1: core of NGC4839 

The origin of this X-ray emission is somewhat unclear. It 
might be the AGN of NGC4839, which is known as a radio 
source, or the emission might be thermal. Fitting a power 
law with the results as shown in Tab.|l| gives a reduced \ 2 
in the order of 1.3. Fitting a thermal plasma improves the 
fit with a reduced \ 2 of about 0.8-0.9. 



5.2. Region 2: the tail of NGC 4839 

This region, for which we only adopted a thermal emission 
model shows high temperature and metallicity. The tem- 
perature is too high for a galaxy, which ranges more about 
1 keV and not around 4.5 keV, as observed. This might be 



region 


center R. A. 


center Dec. 


radius 




(h:min:sec) 


(deg:arcmin:arcsec) 


in arcsec 


1 


12:57:24 


27:29:54 


40 


2 


12:57:17 


27:29:37 


53 




12:57:08 


27:29:06 


65 


3 


12:57:23 


27:27:50 


68 




12:57:15 


27:26:49 


63 


4 


12:56:53 


27:24:27 


205 



Table 2. Definition of the circles for the spectra. The 
coordinates are in J2000. For region 2 and 3 we used two 
circles for extracting the spectra. 



explained by ram pressure stripping: the hot gas is blown 
out as the medium encounters the ICM of the Coma clus- 
ter. Due to the high velocity of the infall of NGC 4839 
into the Coma cluster, shocks could occur, which heat the 
intcrgalactic medium. The fact that we see the tail point- 
ing away from the center of the Coma cluster is a strong 
indication for this scenario (see also discussion). 

5.3. Region 3: the region South of NGC4839 

As can be seen in Fig.^, the region South of NGC4839 
seems to have a different spectrum with respect to center 
or tail of NGC4839. We fit a thermal plasma and a power 
law to the spectrum of this region. The results, again are 
shown in Tab.[l]. Both models for the emission are accept- 
able. The intergalactic medium around NGC4839 could 
be heated by shocks caused from the infall of NGC4839 
onto the Coma cluster (see also the discussion section). 

A power law emission could be linked to the particles 
responsible for the radio emission observed in NGC4839. 
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NGC 4839 is a radio source showing two opposite jets and 
a complex radio tail. This radio tail is extended in S-W 
direction and coincides nicely with region number 3 (see 
Venturi et al. 1990| ). If the X-ray emission in this region is 
linked to the radio tail it could come from Inverse Comp- 
ton scattering of cosmic microwave background photons 
as they interact with the relativistic electrons originating 
from the radio jet. However, in this case it is not clear 
whether XMM has intrinsically enough sensitivity to re- 
solve this IC emission. 

Comparing the surface brightness of region 2 and re- 
gion 3 we find a difference in intensity of a factor of about 
3.5 in the energy band 0.5 — 2.0 keV. If we assume that the 
emission is thermal and thus the radiated energy E cx n 2 
and furthermore that region 2 and 3 are in pressure equi- 
librium we expect a difference in surface brightness of the 
order of 2.5, since region 2 has about kT=5 keV and region 
3 has most likely a temperature of about kT=8 keV. The 
remaining difference between the surface brightness inten- 
sities might, in this case be due to differences in metal- 
licities or extend in the line-of-sight. Region 2 has a high 
metallicity and is in a temperature range, where line emis- 
sion is more dominant than at 8 keV. 

Taking this all together (see also discussion) we think 
that it is more likely that the emission of region 3 is due to 
thermal bremsstrahlung rather than due to IC scattering. 

5.4- Region 4- main body of sub group 

This is the center of the extended X-ray emission of the 
subgroup. The fitted temperature of about kT=4A keV i s 
consistent with ASCA results by Watanabe et al.( 1999| ). 
The metallicity of this group seems to be typical with a 
best fit abundance of 0.3. The abundance seems to be 
significantly smaller than the one observed in the tail of 
NGC4839 (region 2). 



6. Discussion 

6.1. The structure around NGC 4839 

There have been already several studies dealing with the 
physics of galaxies falling onto clusters, like Gaetz et al. 
( |1987| ), Balsara et al. fll994| ) or Stevens et al. ( [1999D (here- 
after SAP99) and references therein. SAP99 found differ- 
ent dynamical features for the infall depending on cluster 
gas temperature, galaxy mass replenishment rate and in- 
fall velocity. SAP99 performed simulations based on a 2D 
hydrodynamic code, which is similar to the one used in 
Balsara et al. ( |1994| ). The ICM follows a beta-model dis- 
tribution with /3 = 2/3 and a mass distribution of the 
which extends up to 32.2 kpc. Similar parameters for the 



galaxy distribution have been used by Gaetz et al. (1987). 
The hardness ratio image Fig.| around NGC 4839 shows 
morphological similarities to their model 2c. The physi- 
cal properties of this model 2c are an intermediate cluster 



at 4 keV and an infall velocity of 1380 km/sec. These 
numbers are slightly different to what is observed in the 
case of NGC4839. Coma has a gas temperature of about 
8 keV and an infall velocity of 1700 km/sec (Colless & 
Dunn 1996, hereafter CD96). Despite the differences in 
actual numbers we interpret the striking resemblances as 
a strong indication for NGC4839 being on its first infall 
onto the Coma cluster (see for example CD96), and not, 



as suggested by Burns et al. (1994) going out of the cluster 



after having already passed through the Coma center. 



6.2. Displacement of NGC4839 and the extended X-ray 
emission 

Another feature which is striking in our observations is 
the displacement of the extended X-ray emission belong- 
ing to the sub group of NGC 4839 and NGC 4839 itself of 
about 8 arcmin (see for example Fig.||). This corresponds 
to roughly 300/i^q 1 kpc. From the optical study by CD96 
NGC 4839 seems to be right in the centroid of the galaxy 
distribution belonging to the sub group. It seems there- 
fore likely that the galaxies belonging to the sub group 
are on average closer to the Coma center than the hot gas 
belonging to the sub group. Such a displacement between 
galaxies and gas is expected when the sub group encoun- 
ters the gas of the main cluster. The ICM of the main clus- 
ter slows down the infall. As the galaxies in the sub group 
show larger densities and therefore a smaller surface-to- 
mass ratio than the gas, the galaxies are less affected by 
the ram pressure force than the intra group gas. In the fol- 
lowing we make a first order calculation to see whether this 
effect is large enough to explain the displacement of about 
300/i^q 1 kpc between NGC4839 and the center of the gas. 
For a rough calculation we estimate that the acceleration 
due to gravitational forces is much stronger than the coun- 
teracting ram pressure force. We furthermore assume a 

(Hughes 

CD96), 



total mass of the Coma cluster of 2 x 10 15 /i 50 1 M 



1989; Bricl et al. 1992 



and that the Coma cluster 
mass is a factor of ten higher than the mass of the infalling 
group (CD96). In a very rough approximation we suppose 
that the Coma mass distribution is to first order point-like. 
This gives an infall velocity v = v / 2GM/r of NGC4839 
of 2400 km/sec, supposing a distance to the Coma clus- 
ter to r — 1.6 M pc (th is is the projected distance see 
also Dow & White |1995[ or CD96 - the most likely infall 



direction is close to the plane-of-sky) . The approximated 
infall velocity is not too different to the calculated value 
of about 1700±^ km/sec by CD96. For the ram pressure 
we use P ra m = pv 2 a n d for the density distribution of 
the ICM a beta-model (p = p (l + r 2 /r 2 )~ 3(3/2 )- For sim- 
plicity we assume (3 = 2/3 and r c — 400h^ kpc, which 
is in good agreement with the ROSAT measurements by 
Briel et al. ( |1992| ). For n e o we use their calculated value of 



3 x 10 _ h 50 cm and transform it into po- For the sur- 
face of the gas encountering the ICM of Coma we suppose 
a circular surface with a radius of 300/i^g 1 kpc (we take 
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the radius of the outermost contour on Fig.||) and a gas 
mass of about 10 13 M©. The force exerted by the ram pres- 
sure is P Tam S = pv 2 S — ma. S is the surface of the sub 
group gas, m its mass and a the deceleration. For the slow- 
ing down we calculate Vdec — J adt — J(PS/mdt/ds)ds 
= S/m J(pVinfaii)ds . We suppose r 2 >> r 2 , which is 
certainly true at distances greater 1.6 Mpc. Filling in all 
parameters we get Vdec = 540 km/sec, which means that 
the actual infall velocity is diminished by 540 km/sec for 
the hot gas of the sub group. Extrapolating a difference in 
velocity between NGC4839 and the gas of 540 km/sec (we 
assume that the ram pressure force is negligible for the 
galaxy) to the time which the group needed to go from 
the virial radius (about 3 Mpc) to the observed distance 
of 1.6 Mpc we see that we obtain a difference in location 
of about 400 kpc assuming a constant infall velocity of 
2000 km/sec for NGC4839. This value which was deter- 
mined with many approximations is of the same order of 
magnitude as the observed displacement between gas and 
galaxy. The differences in location between NGC 4839 and 
the group gas is thus likely to be due to the ram pressure 
exerced by the pressure of Coma's ICM. In addition, more 
detailed modeling of this infall is certainly very interest- 
ing, since it is very likely that the displacement can be 
used to determine the density or temperature of the ICM 
at large radii. 

6.3. The mass ratio NGC'4839 group - Coma cluster 

We obtain a temperature of the sub group of about 
4.5 keV. Comparing this to the central temperature of 
the Coma cluster of about 8-9 keV and applying the 
M-T relation of M oc T 3 / 2 we obtain a mass ratio of 
MNGC4839/M Co ma ~ 0.4. This ratio is very high in com- 
parison to the result of CD96, who find that the sub group 
has 5-10% the mass of the Coma cluster. This discrepancy 
can be either explained by superposition of the gas in the 
group with the hotter ICM of the Coma cluster or the 
fact that the ICM of the NGC4839 group is heated up 
above its initial virial temperature. A plausible scenario 
is that the subgroup is in a process of internal merging. 
- The somewhat disrupted morphology of the extended 
emission seen in Fig.|l| and Fig.|| is an indication for that. 
Such processes of internal merging were already observed 
in hydro dynam ic simulations (see for example Schindler 
& Muller |1993| ). If the ICM of the sub group were heated 
up by a factor of two, which is reasonable, and its initial 
virial temperature was 2 keV, we would obtain a mass ra- 
tio of M nGC mm/M Coma ~ (2/8. 5) 3 / 2 = 0.11, which is in 
good agreement with the values obtained by CD96. 



7. Conclusion 

We performed a first analysis of the XMM-Newton data of 
the NGC4839 group, which has a projected distance of the 
Coma cluster of 1.6/i^q Mpc. In our analysis we resolve 



for the first time a complex temperature structure around 
NGC4839 which can be explained by a first infall of the 
galaxy onto the Coma cluster. 

Furthermore, we measure a displacement between 
NGC4839, the dominant galaxy of the group and the ex- 
tended X-ray emission coming from its ICM. In a simple 
approximation we can explain this displacement by a ram 
pressure force of Coma's ICM acting on the hot gas of the 
sub group as it falls onto the Coma center. 
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